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Abstract 

The BFKL approach for the investigation of semihard processes is plagued by large 
next-to-leading corrections, both in the kernel of the universal BFKL Green’s function 
and in the process-dependent impact factors, as well as by large uncertainties in the 
renormalization scale setting. All that calls for some optimization procedure of the 
perturbative series. In this respect, one of the most common methods is the Brodsky- 
Lepage-Mackenzie (BLM) one, that eliminates the renormalization scale ambiguity by 
absorbing the non-conformal /3o-terms into the running coupling. In this paper, we 
apply BLM scale setting procedure directly to the amplitudes (cross sections) of several 
semihard processes. We show that, due to the presence of /3o-terms in the next-to-leading 
expressions for the impact factors, the optimal renormalization scale is not universal, but 
depends both on the energy and on the type of process in question. 
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1 Introduction 


We discuss the application of the BFKL method jTj to the description of semihard pro¬ 
cesses, i.e. hard processes in the kinematic region where the energy variable s is substantially 
larger than the hard scale Q^, s ^ ^ ^qcd; with Q the typical transverse momentum 

and Aqcd the QCD scale. This approach allows to resum systematically to all orders of 
perturbation series the terms enhanced by leading a”ln"(s/Q^) (LLA) and hrst subleading 
W{s/Q‘^) (NLA) logarithms of the energy. 

In the BFKL approach, relevant physical observables are expressed as a convolution of 
two impact factors with the Green’s function of the BFKL equation. The Green’s function 
is determined through the BFKL equation and is pro cess-independent. The next-to-leading 
order (NLO) kernel of the BFKL equation for singlet color representation in the f-channel 
and forward scattering, relevant for the determination of the NLA total cross section, has 
been achieved in Refs. [ 2 ], after the long program of calculation of the NLO corrections | 3 ] 
(for a review, see Ref. jT]). The other essential ingredient are impact factors, which are not 
universal and must be calculated process by process. Indeed, only a few of them are known 
with NLO accuracy. 

Both the impact factors and the BFKL kernel receive large NLO corrections in the MS 
renormalization scheme. The practical application of this approach to physical processes 
encounters therefore serious difficulties, due not only to the large NLO corrections, but also 
to big renormalization scale setting uncertainties, thus calling for some optimization procedure 
of the QGD perturbative series. 

In this paper we focus on the widely-used Brodsky-Lepage-Mackenzie (BLM) approach | 5 ] 
to face with this problem, which relies on the removal of the renormalization scale ambiguity 
by absorbing the non-conformal /^Q-terms into the running coupling. It is known that after 
BLM scale setting, the QGD perturbative convergence can be greatly improved due to the 
elimination of renormalon terms in the perturbative QGD series. Moreover, with the BLM 
scale setting, the BFKL Pomeron intercept has a weak dependence on the virtuality of the 
Reggeized gluon PE!- 

We apply the BLM scale setting procedure directly to the amplitudes (cross sections) of 
several semihard processes. It is shown that due to the presence of / 9 o"ferms in the NLO 
expressions for the impact factors, the resulting optimal renormalization scale is not universal 
and depends both on the energy and on the type of process in question. We illustrate this 
general conclusion considering in details the following semihard processes: 


• the forward amplitude of production of two light vector mesons in collision of two virtual 
photons, 7*7* —)■ V1V2; 

• the high-energy behavior of the total cross section for highly virtual photons, 7*7* —)■ X; 

• the inclusive production of two forward, high-p^ jets separated by a large interval in 
rapidity Ay (Mueller-Navelet jets), p -|- p —)• jet -|- jet -|- X. 
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At present we do not have a model-independent method to resum the BFKL series beyond 
the NLA logarithms of the energy. Therefore we strictly adhere here to the original formulation 
of the BLM procedure and do not consider its higher-order extensions, such as the sequential 
extended BLM jS] and the Principle of Maximum Conformality [S] (see uni for a review on 
the latter method; see also HHEIEIE] for some recent comparisons between different 
optimization methods). 

The paper is organized as follows: in section 2 we rederive the general expression for the 
NLA BFKL amplitude in the {v, n)-representation; in section 3 we discuss in detail the im¬ 
plementation of the BLM scale setting method, both in exact way and in some approximated 
forms; in section 4 we present the applications of the procedure to the three different processes 
mentioned above; hnally, in section 5 we draw our conclusions and discuss previous studies of 
semihard processes with BLM method. 


2 The BFKL amplitude 


The cross section and many other physical observables are directly related to the forward 
amplitude, which in the BFKL approach can be expressed as follows: 


Im^ (M) = 


(2vr) 


^ /■ (Pq2 

HI J H2 


( 5+200 


^2(—<?2,’So) / 7^ ( — ) G^{qi,q2)- (1) 


( 5—200 


27ri V So 


This expression holds with NLA accuracy. Here, s is the squared center-of-mass energy, 
whereas Sq is an artihcial scale introduced to perform the Mellin transform from the s-space 
to the complex angular momentum plane and cancels in the full expression, up to terms 
beyond the NLA. All momenta entering this expression are dehned on the transverse plane 
and are therefore two-dimensional. $ 1,2 are the NLO impact factors specihc of the process; 
we will see later on three different examples for them. The Green’s function Guj takes care of 
the universal, energy-dependent part of the amplitude. It obeys the BFKL equation 


u}Guj{qi,q2) = - ^ 2 ) + 


d^qK{qi,^G^{q,qi) , 


( 2 ) 


where K{qi,q 2 ) is the BFKL kernel. 

In this section we will derive a general form for the amplitude in the so-called (z/, n)- 
representation, which will provide us with the starting point of our further analysis. We 
will proceed along the same lines of Refs. ng. First of all, it is convenient to work in the 
transverse momentum representation, dehned by 

q\qi) = qM), (3) 

{qM) = S^^\qi - gl) , {A\B) = {A\k){k\B) = j d^kA{k)B{k) . (4) 
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In this representation, the forward amplitude Q takes the very simple form 


Im^ (^) = 




du 


(27r)2 2Tii \so 


$1 


1 


d- - • 


Qi 


<l2 


( 5 ) 


The kernel of the operator K becomes 


^(k,k) = ikikiqi) 

(6) 

and the equation for the Green’s function reads 


i = (a; - k)G^ , 

(7) 

its solution being 

G^ = {uj-k)-G 

(8) 


The kernel is given as an expansion in the strong coupling, 

k = oisk^ + alk^ , (9) 


where 


a.Nr, 


0.0 - 


TT 


( 10 ) 


and Nf. is the number of colors. In Eq. ([^ k^ is the BFKL kernel in the LO, while k^ 
represents the NLO correction. 

To determine the cross section with NLA accuracy we need an approximate solution of 
Eq. dsj). With the required accuracy this solution is 


G^ = {u- + (cn - ) {u - a,K^)-^ + O 


0 ^-l 


-On-1 




( 11 ) 


The basis of eigenfunctions of the LO kernel. 


k^\n,u) = x{n,u)\n,u ), 


77/1 \ / 77 / 1 

X{n,u) = 2 ^lJ{l)-^p ( - + -+iu\-k i- + , (12) 


is given by the following set of functions: 


{q\n,u) = 


TT 


V2 




2 \*^ 2 „in(p 


(13) 


here 0 is the azimuthal angle of the vector q counted from some fixed direction in the transverse 
space, COS0 = qx/Wl- Then, the orthonormality and completeness conditions take the form 


(n', z/) = 


= / ^i{n-n'U _ 


27r2 


= sk — p') Sf, 


(14) 
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and 


i 



du \n, u) (n, u 


(16) 


The action of the full NLO BFKL kernel on these functions may be expressed as follows: 

k\n, v) = v)\n, v) + v) + v) ln(/i^)^ |n, v) 

+ , (16) 

where is the renormalization scale of the QCD coupling; the first term represents the action 
of LO kernel, while the second and the third ones stand for the diagonal and the non-diagonal 
parts of the NLO kernel and we have used 


l3o — 


llA^e 2n/ 


3 3 ’ 

where uj is the number of active quark flavors. 


(17) 


The function calculated in [TB] (see also HZI). is conveniently represented in the 

form 


z/) = (^X^{n, u) - yx(n, z/) - ix'{n, + x(n, z/), (18) 


where 


X{n,u) = -- 


71^ sinh(7rz/) 

2 u cosh^(7rz/) 


(j){n, u) = — dx 


vr^ — 4 


X{n, v) - 6C(3) - x"(n, v) + 2 0(n, v) + 2 0(n, -z/) (19) 


3+1 + ^ 


nf\ 11 + 12 Z /2 


) 16(1 + z/2) 


JnO 


l + 4z/2 

1 + TV7 1 - ^K2 


iV| J 32(1 + z/2) 


X 




1 -I- X 


1 f , f n + 1 


2 IV- 


~ C(2) ) + Li2(x) -|- Li2(—x) (20) 


OO h 


/ \ tJXJ 

+ lnx (V’(n + l)-^/>(l)+ln(l + x) + 5^^^j +5^ 

/ /c=l 


(i-(-i)‘) 


k=0 


- 1 ) 


k-\-l 


k + (n + l)/2 + i^ [V''(*- + n + 1) - V''(*: + 1) + (-1)‘+‘()Z'(*: + n + 1) + ;3'(fc + 1)) 


k + {n + l)/2 + iiy 


— n -M) - ip{k + 1)) 


/?'(-) = I 


v-' 


/ f ^ T 1 




Li2(x) = — dt 


ln(l — t) 
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Here and below = dx{n, i>)/di> and x!'ip-) = dPxip-) v)l,Pv. 

The projection of the impact factors onto the eigenfnnctions of the LO BFKL kernel, he. 
the transfer to the (z/, n)-representation, is done as follows: 


$i(^ 


+ O0 


H-CXO 




— = 'Yh / dv^i{v,n){n,v\qi), 


^2( —^ 2 ) 


+ 00 


+ 00 


^ / ciz/<h 2 (z/,n)(g 2 |n, z/), 


<hi(gi) 1 


$i(z/,n)= I 


<h2(z/,n) = d q2 


Qi 7rV2 

2„ *f*2(-?2) 1 


( 21 ) 


m 


■ 2 „—in4>2 


^2 

The impact factors can be represented as an expansion in Os, 


and 


= o»((Jr) >> 1 , 2 ( 9 ) + aJpR)v'il(q) 
4i,2(«, v) = as(9s) c,, 2 (n, 9 ) + <S,(/ 2 R)c 5 ‘^(n, 12 ) 


( 22 ) 

(23) 


To obtain our representation of the forward amplitude, we need the matrix element of the 
BFKL Green’s function. According to , we have 


(n, z/|G^|n', z/') = 5n,n' 


5{v — v') 


u - a*(/iR)x(n, z/) 
z/) + z/) + f x(^, v) + 2x(n, v) In/i^ + i^))) 


(w - a^(/iR)x(n, z/))2 


(24) 


i^') 


(cj - a^(/iij)x(n, z/))(a; - as{fiR)xin 


,z/0)V^^'^ 7 


Inserting twice the unity operator, written according to the completeness condition (15), 
into (|^, we get 


OO oo 

OO OO f»5+2(DO 


C-AJ n n nQ 


du f s 
2Tri \so 


/ I \ / I I / /\ / / /|^2\ 

X (^|n, z/)(n, z/|Gi^|n,z/)(n,z/|^) , 
9'l '?2 

and, after some algebra and integration by parts, finally 


( 25 ) 


Im* (.4) ^ 


(2r) 


oo 

oo « 

5: jd. 


„ . “sGfllxG.G 

i \ 2 
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«^(AiR)ci(’^, 7 c2(’^, 7 
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(26) 


1 + asifiR) —^ H-^ 

Ci{n,u) C2{n,u) 


/^o 


10 


. d Ci(n, u) 


+a,{fiR)\n— -Ixin^iy) +—x{n,iy) {-xin,i^) + —+ 2\nfij^ + i—\n 

So I oJMc \ C2(n, z/j 


This is our master representation of the NLA BFKL forward amplitude. In the next 
section we will implement on it the BLM scale setting. 


3 BLM scale setting 


The cross section of a process is related, via the optical theorem, to the imaginary part of 
the forward scattering amplitude, 

InioA , , 

' (27) 


a = 


Here we want to discuss the BLM scale setting for the separate contributions to the cross 


section, specified in (26) by different values of n and denoted in the following by Cn- Note 


that the n = 0 case is relevant, e.g., for the total cross sections of 7 * 7 * interactions, Mueller- 
Navelet jet production and the forward differential cross section of the 7 * 7 * —?• V 1 V 2 process. 
Azimuthal angle correlations of produced jets in the Mueller-Navelet process are instead 
associated with non-zero values of n. 

The starting point of our considerations is the expression for Cn in the MS scheme (see 
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Eq. (26)), 


C = 


[2n) 


X 




So 


a;{iiR)ci{n,iy)c2{n, u) 


1 -h asifiR) —^ H-^ 

Ci{n,u) C2{n,u) 


(28) 


+a,{fiR)\n— x(n, v) + 


So 


(do 


10 


. d Ci(? 7 ,, u) 


Q.r-X{n,jy) [-x{n,Jy) + ^+ 2\nn% + ij-\n 

8 Ac V 3 du C2{n,i/) 


In the r.h.s. of this expression we have terms ~ originated from the NLO corrections to the 
impact factors, and terms ~ Q;^ln(s/so) coming from NLO corrections to the BFKL kernel. 
In the latter case, the terms proportional to the QCD /3-function are explicitly shown. For 
our further consideration of the BLM scale setting, similar contributions have to be separated 
also from the NLO impact factors. 

In fact, the contribution to an NLO impact factor that is proportional to (do is universally 
expressed through the LO impact factor. 


n*-^^((f) = v{q) 


/dp 

IK 




(29) 
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where the dots stand for the other terms, not proportional to /do- This statement becomes 
evident if one considers the part of the strong conpling renormalization proportional to n/ and 
related with the contribntions of light qnark flavors. Snch contribntion to the NLO impact 
factor originates only from diagrams with the light qnark loop insertion in the Reggeized 
glnon propagator. The resnlts for snch contribntions can be fonnd, for instance, in Eq. (5.1) 
of [IH]. Tracing there the terms ~ nj and performing the QCD charge renormalization, one 


can indeed confirm (29). 


Transforming (29) to the i/-representation according to (21), we obtain 


~( 1 )/ X /^o 

C. (-.n) = 
c^2\u,n) — 


d / 5 

+*^ci(i/,n) + ( In/r^ + - ) Ci(z/,n) 

d / 5 

+ - I C2{v,n) 


(30) 


and 


^( 1 ) ^( 1 ) 


Cl 


C 2 


/dp 

Wr 




(31) 


It is convenient to introdnce the fnnction / [u), defined throngh 


In ((^ ) =2 [/(n) - In (QiQj)] , 


(32) 


that depends on the given process, where Qi ,2 denote here the hard scales which enter the 
impact factors ci^ 2 ^ The specific form of the fnnction /(i/) depends on the particnlar process. 
According to the properties of the corresponding LO impact factors ( 7 * —)■ R, 7 * —)■ 7 * and 
Mneller-Navelet jet vertex), one can easily check that 


/pp—/•jeti+X+jet2 (^) 0 ) (33) 

for the processes 7 * 7 * —?• X and Mneller-Navelet jet prodnction, whereas for the process 
7 * 7 * (forward electroprodnction of two light vector mesons) this fnnction is not eqnal 

to zero. 


iy) ='ip + 2w) +'ip {3 - 2w) - xIj 




(34) 


Now, we present again onr resnlt for the generic observable showing explicitly all 
contribntions proportional to the QCD /S-fnnction, i.e. also those originating from the impact 
factors: 


C 


n 


1 




al{^iR)ci{n,v)c2{n, v) 


— 00 


'^Here we consider processes whose impact factors are characterized by only one hard scale. This is the 
virtuality of the photon, Q, for the 7 * —> 7 * and 7 * —> impact factors, and the jet transverse momentum, 
|fc|, for the impact factor describing the Mueller-Navelet jet production. 
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1 + asifin) 


c?{n,u) 4"^(n, I/) 00 f 5 


.(1) 


+ 


+ 


ci(n, I/) C2{n,iy) 2Nc V 3 Q 1 Q 2 


—h In 


0R 


+ /H 


(35) 


So 


00 


+a,{nR) In — x{n, y) + 


ANr. 


■—+ - + '°/lr + /(-'' 

2 3 Q 1 Q 2 ' 


where c [^2 = ~ W^e note that the dependence of (35) on the scale fiR is subleading: 

performing in (35) the replacement 




(36) 


one indeed obtains the same expression as before with the new scale at the place of the 
old one fiR, plus some additional contributions which are beyond the NLA accuracy. 


As the next step, we perform a finite renormalization from the MS to the physical MOM 
scheme, that means: 


^MS ^ ^MOM 


1 + 


a 


MOM 


TT 




with T = 


{1 + -I] , 


(37) 


(38) 


rj^coni 00 A 




where I = —2 ~ 2.3439 and ^ is a gauge parameter, fixed at zero in the following. 


Inserting (37) into (35) and expanding the result, we obtain, within NLA accuracy, 

afO'^{fj.n)x{n,u) 


c, 


MOM 


ci(n, z/)c2(n, v) 


00 


X 


l+df°^(/i«) 


0^\n,v) C2^^(n, I/) 

ci(n, I/) C2(n,i/) Nc 




2Nr V3 


Q 1 Q 2 


+ («™(/^r))^ In ^ <1 X(n, v) + ^^X(n, v) 


“'conf 


AN, 


x{n,y) 


0R 


Xin.p) , 5 , 


N,. 


+ fi^) “ ^ 1 ^ 3 -^ 


(39) 


The optimal scale is the value of fiR that makes the expression proportional to 0o vanish. 
We thus have 


Cf = 


(27r) 


du ( — 
So 




(^M0M(^BLM))3 


— 00 

































X Clin, z/)c 2 (n, v) 


/^o 

Wr 


,BLM\2 


3 Q 1 O 2 ' 


211 + -/ 


(40) 






(/^ 


iJ 


So 


BLM\2 


^+/(.)-2(l+?/ 


= 0 


In the r.h.s. of (40) we have two groups of contributions. The hrst one originates from the 


/9o-dependent part of NLO impact factor (29) and also from the expansion of the common 
cx^ pre-factor in (35) after expressing it in terms of The other group are the terms 

proportional to d^^lns/sp- These contributions are those /3o-dependent terms that are 


proportional to In s/sq in (35) and also the one coming from the expansion of the 
factor in (35) after expressing it in terms of cx^ 


MOM 


The solution of Eq. (40) gives us the value of BLM scale. Note that this solution depends 


on the energy (on the ratio s/sq). Such scale setting procedure is a direct application of 
the original BLM approach to semihard processes. Finally, our expression for the observable 
reads 


^BLM _ 


(2ir) 


du I — 
So 


(41) 


X («f°^(hfl^“))^ci(n,i/)c2(n,z/) 




M), 


n, u) 




+ 


Ci(n, z/) C 2 (n, z/) 


^ + 


Nr 


where we put at the exponent the terms Ins/so, which is allowed within the NLA 

accuracy. 


Unfortunately, Eq. (40) can be solved only numerically, thus making the scale setting a 


bit unpractical. For this reason, we will work out also some analytic approximate approaches 
to the BLM scale setting, which have the merit of a straightforward and simple application. 
We consider the BLM scale as a function of u and chose it in order to make vanish either the 
hrst or the second (~ d^^^lns/so) group of terms in the Eq. (40). We thus have two cases: 


case (a) 


^BLM,a _ 


= ' 5 iQ 2 exp 


2|i + h)-/(0-5 


(42) 


(2^)= 


dv I — 

So 




X Ciin,u)c 2 in,u) 


X 


1 +<; 

I Ci(n, z/) C 2 (n, z/) 


c^i \n,h') C 2 ^^(n, z/) 


+ 


Nr 


(43) 


case (6) 


(/i 


= QiQ 2 exp 


2(1 + -/ 


/(^) - ^ + \xi^,n) 


(44) 
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BLM,b 


(2vr) 


dv ( - 


X 


X -j- cn 


MOM/ .BLM 


{f^R,b 


X (af°^(/^S'^))^Ci(n,z/)c2(n,i/) 

c^^\n, v) ^ 4'^(n, v) ^ ^ —x(n, z/) 


(45) 


ci(n, z/) C 2 (n, z/) 


iV. 




The other possible option for the BLM scale setting could be related with the requirement 


that the entire expression in the integrand of (40) vanishes, which leads to the following 
condition 


case c 


5 , 

—h In 
3 


(/^ 


/?,C 


Q 1 Q 2 


+ /(^) ~ 2 




^,MOM(^BLM^ In 


a 




so 


1 "1“ 


MOM(/,BLM)i^^ 


-R,c 


X{'n,v) 


(46) 


One should mention, however, that such approach to the BLM scale setting has a limited 
applicability, since the denominator in the r.h.s. of ( [4^ vanishes at some value of z/ = z/, 
given by 

_ <? A/I r) f/I 

(47) 


l+af°"ln/h!;Tl = 0, 


So 2 

which prevents us from dehning /i^^^(z/) in the entire z/ range. Nevertheless, one can try to 
use such method in those cases when the product of the two LO impact factors ci(rz,, z/)c 2 (n, v) 


is a function decreasing so rapidly to guarantee the convergence of the z/-integration in (41) 


in the z/-region where there is no problem with the solution of Eq. (46) 


Note also that all three approaches to BLM scale hxing discussed above, and given in 


Eqs. (42), (44) and (46), could be applicable only to processes characterized by a real-valued 


function /(z/). For some processes this is not the case. In particular, the inclusive production 
of two identihed hadrons separated by large interval of rapidity in proton-proton collisions, 
p + p —)■ hi + h 2 + X, is described by a complex-valued function, /*(z/) = /[—v). This 
can be easily seen calculating f{u) from Eq. (77) of [12] for the identihed hadron production 
impact factor. In such cases one can use only the BLM scale hxing method which relies on 


the numerical solution of Eq. (40). 


4 Applications 


In this section we apply the BLM approach to a selection of semihard processes. For the 
energy variables we will use notations 


Y 



Yo = 



(48) 


In our numerics we use the following settings: Hf = 5 and as{Mz) = 0.11707 for the number 
of active havors and the value of the strong coupling. 
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Figure 1 : Left: BLM scales for the process 7 * 7 * — )• V1V2 (see the text for details). Right: Forward 
amplitude for 7 * 7 * —)> V1V2 at Iq = 0 . 


4.1 Electroproduction of two vector mesons 


We start with the description of the forward amplitude for the production of a pair of light 
vector mesons in the collision of two virtual photons, 7 * 7 * —^ViV 2 . Such processes could be 
studied in experiments at future high-energy e’''e“ colliders, see [SniEQEg for estimates of 
the cross section in the Born approximation. The BFKL resummation for these processes was 
considered in |23], where the inclusion of NLO effects was limited to the corrections to the 
BFKL kernel. In the papers na, some of us performed a complete NLA BFKL analysis for 
the forward amplitude of these processes, including the NLO corrections also to the 7 * —>■ K 
impact factors |24j . Very large NLA corrections to the forward amplitude were found, therefore 
in ng the Principle of Minimal Sensitivity (PMS) j25] approach was used to optimize the 
perturbative series. 

Here we present numerical results for the forward amplitude obtained with the BLM 
optimization method described above. 


We consider the case of equal values of photons virtualities, Qi = Q 2 = Q, and, following 
the first of Refs. ng, present our numerical predictions for the forward amplitude multi¬ 
plied by some kinematic factors, Im^ (M) Q‘^/{SD 1 D 2 ), calculated at Q = 50 GeV, where the 
expressions for D 12 are given in Eq. (14) of the first of Refs. |15| . 


For the considered process only the n = 0 term contributes and the f{u) function is given 
in (34). We will try all approaches to the BLM scale setting described in previous section. 
In particular, for this process the product of two LO impact factors ci{n, z/)c 2 (n, u) vanishes 
very fast for \v\ > 1 , therefore in the relevant integration z/-range, \v\ < 1 , we can find the 
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solution of Eq. (46) and determine the BLM scale as a function of v and energy. 


In Fig. Bi eft) we show the values of the BLM to kinematic scale ratios, as 

functions of E — Yq, obtained in four different cases. By “exact” case we denote the scale 
obtained solving numerically Eq. (40) for each value of ln(s/so) = b' — Iq- In the other 
three approaches, the BLM scales depend on u: the scales for cases (a) and ( 6 ) are given by 
Eqs. ( |4^ and (44), respectively; the case (c) corresponds to the numerical solution of Eq. (46) 
for each values of u and Y — Yq. The i/-dependent scales, cases (a), ( 6 ) and (c) are shown in 
Fig Bl eft) for the particular value of i/ = 0 . 

Approximate approaches to the scale setting give energy-independent BLM scales (see 
cases (a) and (b) in Fig. [^left)), whereas an exact implementation of the BLM rule leads in 
general to the scales which depend on the energy of the process (see cases (c) and “exact” 
in Fig. I^left)). In fact, the approaches (a) and ( 6 ) can be considered as a low- and a high- 
energy approximation to the case (c), where the BLM scale setting prescription is implemented 
precisely. 


Nevertheless, as we already mentioned above, the condition (46) could not be resolved 
for all processes. Therefore we defined also a method which could be universally applied 
and which we call here “exact”. It gives a i/-independent BLM scale and it is based on the 


requirement of vanishing of the integral in Eq. (40), contrary to the approach (c), where we 
try to make vanish the integrand of the same equation for each separate value of z/. 

In Fig. [fright) we show our predictions as functions of the energy for the forward ampli¬ 
tude calculated with all the four different methods described above: cases (a) and {b) were 


calculated using Eqs. (43) and (45), cases (c) and “exact” using Eq. (41) with the correspond¬ 


ing choices of the scales. The result of the BFKL resummation depends not only on the 
renormalization scale fiR, which is fixed here with the BLM method, but also on the energy 
scale So or Yq. In Fig. [fright) we present the results obtained with the choice of this scale 
dictated by the kinematic of the process, Sq = or Eq = 0 - A more reliable estimation could 
result from fixing the value of Yq according to some optimization method, such as PMS, but 
this goes beyond the scope of present paper. 

As we can see in Fig. [fright), our predictions obtained with precise implementations of 
BLM method he inbetween those derived with the use of the two approximate realizations. 
Note that the difference between the two explicit methods, cases (c) and “exact”, is sizeable 
and increases with the energy. This is related to the fact that these two approaches are not 
equivalent, and the scales in the case (c) are larger than those in the “exact” one. Note also 


that, with the growth of energy, the value of z/ = z/ where the solution of Eq. (46) has a sin¬ 


gularity decreases, see Eq. (47), and approaches the z/-range important for the determination 
of our observable. 


4.2 7 * 7 * total cross section 

In |2n] some of us studied the 7 * 7 * total cross section in the NLA BFKL approach consid¬ 
ering two different optimization methods of the perturbative series. One of them was the BLM 
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Figure 2 : BLM scales for the process 7*7* —X (see the text for details). 


method, cases (a) and (b), described above, where the Eqs. (43) and (45) were transformed 
back to MS scheme. In that paper we fixed the photon virtualities and, correspondingly, the 
number of active flavors nf in order to make a comparison with LEP2 experimental data. 
Here, we are interested in the general features of the BLM scale setting procedure, therefore 
we prefer to fix the photon virtualities as in the two-meson production; Qi = Q 2 = Q = 50 
GeV with n/ = 5. 

In Fig. as in the case of the vector mesons, we show the four different ratios 
versus V — Yq. The four cases (a), (b), (c) and “exact” are defined exactly as in the previous 
subsection. As we mentioned before, cases (a) and (b) are independent on the energy of the 
process, but depend on the kind of process through the f{u) function. In particular, for the 


production of a pair of light vector mesons the function is given by Eq. (34), while for this 


process it is /{u) =0 (see Eq. (33)). 


For this process we only discuss here the BLM scale setting and do not present its cross 
section. The 7 * 7 * cross section was already considered in |2S], where serious problems were 
found, related with the very large values of NLO corrections j2Z] for the virtual photon impact 
factor. For details and an extended discussion of this issue, we refer the reader to (26]. 


4.3 Mueller-Navelet jets 

The last semihard process that we consider is the production of two forward high-pj’ jets 
produced with a large separation in rapidity Ay (Mueller-Navelet jets [2H])- Such process 
was studied at Large Hadron Collider (LHC): the CMS collaboration provided with data for 
azimuthal decorrelations ra that can be expressed, from a theoretical point of view, by ratios 
Cm/Cn, where Cn are to be averaged over kj. (jet transverse momentum) and yJ^ (rapidity jet). 
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Then, in order to match the kinematic cuts used by the CMS collaboration, we have 

^yi,max ^y2,inax pOO POO 

Cn= dyi / dy 2 / dkj^ / dkj^6 {yi - y 2 -Y)Cniyj^,yj 2 ,kj^,kj^) ,(49) 

'^yi,min ^ y2,niin Jkj^ 

,min Jkj^ 

,inin 

with 2 / 2 ,min 4.7, l/l^max 2/2,max 4.7 Q aild fcj]^,min ^J 2 ,min 35 GeV. The 

comparison between experimental results for jets with cone radius R = 0.5 produced at a 
center-of-mass energy of \/s = 7 TeV and theoretical calculations was done in |3l], where the 
exact NLO impact factors calculated in |32) were used, and in |30|, where the NLO impact 
factors were taken in the small-cone approximation as calculated in j33] 


In this section we use the same kinematic settings as in [30] and present the BLM 
scale setting for the Mueller-Navelet jet production. In particular, we consider the ratios 
i^T'/Vk7, kj^ as functions of T — ho for ri=0, 1, 2 and 3 and recall that, for this process, 
the function f{u) = 0 is zero. The results are shown in Fig. were the three lines, violet, 
green and blue, denote the cases (a), (6) and “exact”, respectively. For this process it is not 
possible to consider the case (c) because the product of LO impact factors ci(n, i/)c 2 (n, v) is 
a function that does not decreases so fast, so that the i/-interval needed for the integration 
includes the value V, defined by Eq. (47), where the method is not applicable. 


Due to the integration over the jet variables kj^^^ and fhe derivation of the “exact” 
curve here is a little bit different from that of the other two processes. In this case, in order 
to get the ratios /^blm/ ^kj^kj^, we write = nriR^Jkj^kj^ and we look for rriR such that 
Eq. (40) is satisfied. 


On the contrary, since cases (a) and (6) are independent of the energy of the process for 
n = 0, the two curves, “BLMa” and “BLMb” are equivalent to those in Fig. since also in the 
present case /(z/) = 0. Moreover, note that, for n = 1, x(n = 1, z/ = 0) = 0 and therefore in 
Fig. I^the curve “BLMb” overlaps exactly the curve “BLMa”. 


In Fig. 1^ we present some ratios Cm/Cn versus Y, where we make use of the scales shown in 
Fig-i In all cases shown in Fig. the factorization scale fip entering the MSTW2008nlo [35] 
parton distribution functions was chosen equal to the renormalization scale /i/j and the BFKL 
energy scale Fq was fixed at zero. One could look for optimal choices of the scale Fq, based 
on the PMS method, for instance, but this goes beyond the scope of the present paper. 
The results shown in m are a little bit different from those shown here, because there we 
transferred back all formulas to the MS scheme. Moreover note that now we have an extra¬ 
curve (BLMexact) fo which fipi was obtained solving Eq. (40). 


5 Summary 


In this paper we have focused on the BLM method to set the renormalization scale in a ge¬ 
neric semihard process, as described in the NLA BFKL approach in the {u, rz,)-representation. 

^In |30| it was mistakenly written yi^min = 0, although all numerical results presented there were obtained 
using the correct value = —4.7. 

^ For a critical comparison of the different expressions for the forward jet vertex, we refer to |34| . 


14 









We found that the BLM scale setting procedure is well defined in the context of semihard 
processes described by the BFKL approach within NLA accuracy. The straightforward appli¬ 


cation of the BLM procedure leads to a condition to be fulfilled, Eq. (40), which defines the 


optimal renormalization scale depending on the specific process and on its energy. Our main 
observation here is that, due to the presence of /3o-terms in the next-to-leading expressions 
for the process dependent impact factors, the optimal renormalization scale is not universal, 
but turns out to depend both on the energy and on the type of process in question. The 
non-universality of the BLM scale setting in exclusive processes was observed already in 


Note that the above-mentioned ~ /3o-contributions to NLA impact factors are universally 


expressed in terms of the LO impact factors of the considered process, see our Eqs. (29) 


and (30). Thus, they could be easily calculated for all processes, even in the case when the full 


expressions for the NLO corrections to the impact factors are not known. Such contributions 
must be taken into account in the implementation of BLM method to the description of cross 
sections of semihard processes, because all contributions to the cross section that are Po 
must vanish at the BLM scale. 


Such an “exact” implementation of the BLM method could be difficult since it calls for 
the solution of an integral equation, Eq. (40), for each value of the energy of the process. 
This equation can be solved, in general, only in numerical way. Therefore, we considered 
several approximated approaches to the BLM scale setting. One of them, the closest to the 
“exact” one and labeled (c), consists in imposing the vanishing of the integrand appearing in 
the above-mentioned general condition and leads to an optimal BLM scale depending also on 
the //-variable. This approximated method has a validity domain in the //-space and can be 
applied only if the relevant range of the //-integration giving a physical observable falls inside 
this validity domain. Other approximated approaches, labeled (a) and (6), can be viewed as a 
sort of low- and high-energy approximation of the case (c) and of the “exact” determination. 


We have compared these different approaches in the study of the total cross section and 
of other physical observables related with the forward amplitude in processes such as the 
electroproduction of two light vector mesons, the total cross section of two virtual photons and 
the production of Mueller-Navelet jets|^ Note that the formulas for the approximate cases (a) 
and (6) were already used by us, without derivation, in our recent papers [261 l30] . Here 
we presented in full details the implementation of the BLM method for arbitrary semihard 
processes, considering both its exact and approximate forms. 


We could observe that, in general, the BLM scale setting in the cases (a) and {b) provides 
with a range inside which lie the “exact” the case (c) determinations. This is not the case for 
the Mueller-Navelet jet production where, as discussed in the text, due to some peculiarities 
in the definition of the observables imposed by the experimental cuts, the natural ordering 
between the optimal scales in the cases (a), (6) and “exact” is sometimes lost. It turns out, 
however, that azimuthal correlations and ratios between them in the Mueller-Navelet case are 
less sensitive to the different approaches to BLM scale setting than in the other two processes 
considered in this work. 

'*^For all these cases the expression of the amplitude is known within the NLA as convolution of NLO impact 
factors with the NLA BFKL Green’s function. 
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Note that previous applications of the BLM method to the description of 7 * 7 * total cross 
sections |371 [381 m [391 relied on the use of LO expressions for the photon impact factors. 
In [371 [38] the 7 * 7 * total cross section was considered in LLA BFKL, since the NLO corrections 
to the BFKL kernel were not yet known. However, in ISZllMI the /do-part of the first correction 
to the Born amplitude (he. the t-channel two-gluon exchange) was considered in order to 
establish the renormalization scale. Such approach to the scale setting is closely related to 
our case (a) (scale fixed from the correction to the impact factor). Indeed, considering the 
expansion of the BFKL amplitude (26), one can see that the first, ~ correction to the Born 
amplitudes originates entirely from NLO parts of the impact factors. Comparing Eq. (5.5) 
in |3H| with our Eq. (42) for f{i>) = 0 , as appropriate for the 7 * 7 * process, one can see that 
they agree except for the term that, in our approach, derived from the change to the MOM 
scheme. One can therefore refer to [3Z[[38| as to the first (approximate) application of the 
BLM scale setting to a BFKL calculation. In [7[ [39] the 7 * 7 * total cross was considered using 
full NLA BFKL kernel, but with LO approximation for the photon impact factor. In respect 
with the BLM scale setting, such approach is equivalent to our approximate case {b). 


In [3l] the BLM method was applied to Mueller-Navelet jet production; although the 
full NLO expression for the jet impact factor was used, the above-discussed effect of /do- 
contributions to NLO jet impact factors on the choice of the BLM scale was overlooked. 
Therefore in |3T] the value of the BLM scale which was obtained is similar to the one used 
in I3 W\ and, as such, coincides with our approximate case ( 6 ). Our results presented in 
Fig. a allow to assess the inaccuracy in BLM predictions for different Mueller-Navelet jet 
observables related with approximated approaches to the BLM scale setting. 


In conclusion, the BLM method for scale setting, which was proposed more than three 
decades ago on a strong physical basis, remains a fundamental tool for perturbative cal¬ 
culations and has lead to many successful comparisons between theoretical predictions and 
experimental data. In this paper we have provided with the general paradigm for its system¬ 
atic application to an important class of processes, i.e. semihard processes within the NLA 
BFKL approach, thus filling some gaps left open by previous approximated or incomplete 
approaches. We believe that this will increase the future significance of the method. 
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Figure 3: BLM scales for Mueller-Navelet jets (see the text for details). 
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Figure 4: Azimuthal decorrelations for Mueller-Navelet jets (see the text for details). 
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